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Abstract. This paper is based on recent work which provided an exact analytical 
description of scattering fidelity experiments with a microwave cavity under the 
variation of an antenna coupling [Kobcr et al., Phys. Rev. E 82, 036207 (2010)]. 
It is shown that this description can also be used to predict the decay of the fidelity 
amplitude for arbitrary Hcrmitian perturbations of a closed system. Two applications 
are presented: First, the known result for global perturbations is re-derived, and 
second, the exact analytical expression for the perturbation due to a moving S-wave 
scatterer is worked out. The latter is compared to measured data from microwave 
experiments, which have been reported some time ago. Finally, we generalize an 
important relation between fidelity decay and parametric level correlations to arbitrary 
perturbations. 
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1. Introduction 

During the last decade, approximately, considerable efforts have been dedicated to 
the quantitative prediction of the fidelity decay in chaotic /diffusive quantum systems 
and classical wave systems [TJ EJ [3] H] (see also j5] and references therein). A very 
successful approach has been based on random matrix theory, adopting the so called 
Bohigas-Giannoni-Schmit conjecture [6]. Applied to the current setting, it suggests 
that quantum systems with chaotic classical counterpart ("chaotic quantum systems" 
for short) as well as diffusive wave systems show a universal response to perturbations 
which can be calculated within an appropriate random matrix model [7] . The first exact 
analytical results in this respect have been obtained by Stockmann and Schafer [8j [9] 
using super-symmetry techniques similar to those for the calculation of correlation 
functions between scattering matrix elements in [TO]. More recently, exact analytical 
results have also been found for scattering systems, where the fidelity amplitude, an 
expectation value, is replaced by the "scattering fidelity", which is a product of two 
transition amplitudes [III [12]. These results, published in [13], have been obtained by a 
simple but powerful modification of the Verbaarschot-Weidenmiiller-Zirnbauer ("VWZ" 
for short) formula from [TO] . We call this approach the "scattering approach to fidelity". 

As shown first by Kohler et al. [T3], for a global perturbation of a completely 
diffusive system, the fidelity amplitude can also be calculated from the parametric 
level correlations. Subsequent generalizations have been discussed in [15], [16] and [17] . 
Originally, parametric level correlations have been introduced in the area of disordered 
systems with diffusive dynamics [18] [19]. At that moment, they have been considered a 
universal signature of chaotic /diffusive dynamics, where the functional does not depend 
on the perturbation applied. However, in [20J it was shown that certain types of 
perturbations may lead to pronounced deviations. Thereby, it became clear that the 
"universal" theoretical prediction of Simons and Altshuler [19] only applies to global 
perturbations, not local ones, where the perturbation operator has only a few eigenstates 
with non-zero eigenvalues. The perturbation due to the displacement of a small scatterer 
discussed in [20] is precisely of that latter type. 

The purpose of the present paper is twofold. Firstly, we use the scattering 
fidelity approach from [13] to derive exact analytical expressions for the fidelity decay 
of chaotic/diffusive wave systems in the presence of completely general Hermitian 
perturbations. This allows us, to re-derive the known result for the decay of the fidelity 
amplitude due to a global perturbation [8] [9]. In the second example, we use it to 
derive an exact analytical expression for the decay of the fidelity amplitude due to the 
displacement of a S-wave scatterer (local perturbation). In [21] . the decay of the fidelity 
amplitude has been obtained from experimental for such a case [20J . 

Secondly, we generalize the relation between the fidelity amplitude and the 
parametric level correlations from [T3], to arbitrary perturbations. To do so, we compare 
the analytical expression for the parametric level correlations [22] and its analogue for 
the fidelity amplitude for general perturbations. Our result is important, as it allows 
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to calculate the fidelity amplitude for spectral data, only. It thereby shows that the 
fidelity amplitude is a basis independent - which is surprising taking into account that 
the perturbation may be completely arbitrary. From a practical point of view, one 
may easily find situations, where the measurement of level spectra and their variation 
under certain perturbations is easier and more accurate than any fidelity measurement. 
The new relation shows, that a measurement of parametric level correlations provides 
exactly the same information about an a priori unknown perturbation than a fidelity 
measurement. 

The present paper is organized as follows: In the next section, we follow [13] to 
describe the connection between scattering fidelity [H] and scattering matrix correlation 
functions as considered in [10]. We use this connection to derive an exact analytical 
expression for the fidelity amplitude valid for arbitrary perturbations. In section |3] we 
discuss the differences between local and global perturbations, and we use our general 
formula to re-derive the known result for a global perturbation. In the remaining part 
of that section, we calculate the fidelity amplitude in the case of a moving scatterer 
and compare the resulting theoretical prediction to experimental data from [2T). In 
section HJ we evaluate the general integral expression for the fidelity amplitude in the 
perturbative/long time limit. In section |5]we generalize the relation between parametric 
level correlations and the fidelity amplitude to arbitrary perturbations. Conclusions are 
presented in section El 

2. Scattering approach to fidelity 

In this section, we introduce the central quantity of this work, the fidelity amplitude of 
a closed quantum or classical wave system, with quantum chaotic or diffusive dynamics. 
We assume that random matrix theory can be used to describe the fidelity decay. While 
the first part contains some general statements about fidelity and the random matrix 
models used, the second part describes the description of the algebraic scattering model 
to which the fidelity problem is mapped. This mapping, introduced in [13] provides an 
exact analytical description of the fidelity decay. 

2.1. Fidelity 

The fidelity and the fidelity amplitude for a Hamiltonian H a = H + W a with 
perturbation W a are defined as 

F(t) = \f(t)\ 2 , f(t) = (a\e 2niH ^e- 27riH ^\a) (1) 

where \a) is the initial state and W a is the perturbation depending on an external 
parameter a. We assume that the energy is measured in units of the mean level spacing 
do in the spectrum of H , and time in units of the Heisenberg time in = "2TTh/d Q . As a 
result, the variable t in ([T]) becomes dimensionless. 
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For our purpose it will prove convenient to write the perturbation in terms of a 
normalized eigenbasis: 

Wa = X>» \vMvjl , (2) 
3 

where the orthonormal eigenstates { } are assumed to be independent of a. This 
is normally well fulfilled in the case of global perturbations, and also in the case of 
many types of local perturbations, such as point like scatterer. A detailed discussion 
is given in section [3j In other words, (jSJ) implies that [W a , W@] = for any a, (3 in the 
allowed range. Note that it is often possible to consider H + W a as the unperturbed 
Hamiltonian H' Q . Then, W' a = and Wp = Wp - W a such that \W' a , W' p ] = 0, trivially. 

Returing to our original setup, we choose Hq from one of the invariant ensembles, 
the Gaussian orthogonal ensemble (GOE) or the Gaussian unitary ensemble (GUE) [23J. 
Correspondingly, we assume that the perturbation W a can be diagonalized either by an 
orthogonal (GOE case) or unitary (GUE case) transformation. In either case, we arrive 
at 

H a = H + Y^Wj(a) , (3) 

3 

without changes in the random matrix ensemble for H . Here, the states \i) simply are 
the elements of the canonical basis of a complex vector space C^, where N may be 
assumed arbitrarily large but finite. In this situation, one may use the results of [13] to 
calculate the fidelity amplitude averaged over H as the average of a certain scattering 
matrix correlation function within the framework of statistical scattering [10]. In what 
follows, we concentrate on the GOE case. The GUE case (which turns out to be even 
simpler) may be treated along similar lines, using [21] . 

2.2. Scattering matrix correlation functions 



According to [TO], the scattering matrix may be written as 

Sab(E) = S ab - 2m Ut 1 V , (4) 

— -"cff 

where H e g = H — in VV^ with H from the GOE. The indices a, b denote the 
scattering channels which may be chosen in such a way that the column vectors of the 
rectangular matrix V are orthogonal. Analogous to ([3]) it is thus possible to diagonalize 
the perturbation such that 

H cS = H - in |a)(a| (5) 

a 

where the parameters 7 a > are the real and positive eigenvalues of V (eigenstates 
corresponding to zero eigenvalues are ignored) . According to [10] , the average S-matrix 
(averaged over Hq), is then given as 



E(S ab ) = , n a = ^ , (6) 
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where we have assumed that the average level spacing for H is equal to one. Here, 
we introduced the somewhat unusual notation E( . . . ) for the ensemble average over 
the Gaussian random matrix ensembles, to avoid possible conflicts with the Bra-Ket 
notation used below. The main result of [10] consists in a triple integral which gives the 
spectral correlation function between different S-matrix elements 

C[S* ab , S cd }(w) = E[ Sab(E)* S cd (E + wd)}- E[S ab (E))*} E[S cd (E + wd)} , (7) 

depending on the transmission coefficients T a = 4n a (1 + n a )~ 2 i only. Due to the 
convolution theorem, the Fourier transform of these correlation functions yields an 
average over different amplitudes of the evolution operator for the effective Hamiltonian 
if e ff. Namely, for t > 0: 

C[S* ab , S cd }(t) oc (S ab (t)* S cd (t)) = E( (6|e 2 ^ ff *|a)(c|e- 2 ^ *|d) ) , (8) 

It is still assumed that the Hamiltonian is written in the eigenbasis of the coupling VV'. 
Therefore, the states \a), \ b), \c) and \d) represent elements of the canonical basis in C . 
In what follows, we will only be concerned with the case when c = a and d = b. This 
yields for the correlation function in (jSJ): 

C[S* ab , S ab }(t) = 6 ab T 2 (Z J 2 )j + (1 + 5 ab ) T a T b (ZP ab )j , (9) 

where the angular brackets (. . .)/ denote the following weighted double-integral: 



max(0,t— 1) 



(2u + l)(t 2 -r 2 + x 2 



and where we have introduced the following short hands: 

nl-TJt-r) , 9 2r + l 

A — , x 2 = u 2 — — . (11) 

j= i Jl + 2T 3 r + T 2 x 2 2u + 1 



Note that changing the integration variable from u to x yields 
2u 2u 2 



2xdx = (2r + 1) 



du dx , . 

du => = ; , 12 

2u + 1 Vx 2 + 2r + 1 



_2u + l (2u + 1) 2 
such that ( ITOl) may be written equivalently as 

./m.x(0,,-i) Jo Vx 2 + 2r + 1 (t 2 - + x'f 
That expression can be compared directly to the results in [HI [9]. 



2.3. Connection to fidelity 

Starting from (CQ), we insert the projection onto a random state \b)(b\ into the definition 
of the fidelity amplitude: 

/(*) f ab (t) oc E( (ale 2 ^ 1 \b)(b\ e" 2 ^* |o> ) , (14) 

where the average over the random state \b)(b\ simply yields the identity times a 
normalization constant equal to the inverse Hilbert space dimension. As a result, 



Scattering approach to fidelity decay 



6 



we obtain the product of two transition amplitudes which may be considered as a 
scattering fidelity as introduced in [TT] . Comparing the effective Hamiltonian (jSJ) with 
the perturbed Hamiltonian for a closed system as given in we find that they share 
the same structure, and that we only need to allow the coupling parameters 7 a to become 
complex to unify both descriptions. 

In [13], it has then been shown, that the analytical result for the correlation 
functions in (Q can be generalized to different effective Hamiltonians H e g and H' cS , 
which differ only in the eigenvalues 7 a and ■y' a . In that case, one just needs to calculate 
the effective transmission coefficients 

2 «• + K*) 

Ta = - ^ / w 3 ' r 15 

3 (1 + ^(1 + *;) 1 ; 

from the coupling parameters n a (corresponding to H e g) and n' a (corresponding to 
H' eS ), as defined in (JBJ). Then, the double integral in (|TU|) yields the scattering fidelity, 
when replacing the transmission coefficients in the term Z by the effective transmission 
coefficients defined in (fl5|) . 

Restricting ourselves to closed systems with Hermitian perturbations, we obtain 
from the comparison of ([3]) with flH]) that — VK^fj = Wj(a), such that according to ([6]) 



7T 



7^ i Ti Wj (a) , in Wj (j3) 



Finally, to make sure that we really have a closed systemxs, we need the transmission 
coefficients T a and TJ, to be negligibly small. This means that the dynamics of the system 
is probed from the outside via scattering channels which are so weakly coupled to the 
system, that their effect on the dynamics is negligible. The functions to be integrated 
in 09]) then become 

J a ->2t 

P ab ^ P = 2[r 2 + (2r + l)t-t 2 -x 2 ] . (17) 

Thereby, we obtain for the scattering fidelity 

fMKj}, {«;} ; t) oc 6 ab T 2 At 2 (Z)j + (1 + 5 ab ) T a T b (Z P >/ • (18) 

Here, we indicate explicitly the dependence of the scattering fidelity on the coupling 
parameters {ftj} and {ft'j} as their value will become important below, where we discuss 
normalization. 



Normalization In order to calculate the fidelity amplitude from the scattering fidelity 
fab{t), there is still the problem of normalization to be solved. This is because f ab (t) 
becomes an auto correlation function for H c q = H' cS , which still decays to zero in time, if 
the coupling to decay channels is finite. In [TT] . this problem has been solved by dividing 
the scattering fidelity through the geometric mean of the auto correlation functions of 
H c {[ and H' eS . Below, we will see that this normalization procedure is somewhat simpler 
in the case of closed systems. 
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As mentioned earlier, when one is really interested in fidelity and \b)(b\ has been 
inserted just for convenience as described in (JUJ), one can normally assume that a ^ b. 
In addition, the case a ^b arises in the case of an explicit scattering fidelity experiment, 
where in- and out-going channels are chosen to be different (transmission measurement). 
Then, the formula for f a b{t) simplifies to 

f ab (t) cx T a T b (P Z)i . (19) 



In order to apply the normalization scheme from [IT] , we note that for the auto 
correlation functions: 

/«*({«;}, {«;} ; t) = / a6 ({4>, {kJ} ; t) cx T a T b (P )j . (20) 

This follows from the fact that Ka + k* = k'- + «;'•* = since the coupling parameters are 
purely imaginary in both cases. That implies that the effective transmission coefficients 
are zero, so that Z becomes equal to one. Since the auto correlation functions are the 
same, the geometric mean is also the same, and 

t (f K i (VI ■ +\ - TaTb ( P ° Z ^ - ( P ° Z ) 1 con 

l a ±b \Pq)i \Pq)i 

Now, it has been shown in [8] that for Z = 1, the resulting double integral yields 

7m«(0.,-i) Jo Vx' + 2r + 1 (A - r' + x 2 ) 2 

for any t > 0, so that 

/a6({s} ; {4}^) = ( P 0^)7- (23) 

This formula constitutes the first important result of our work, since it gives an exact 
analytical expression for the fidelity amplitude of a chaotic /diffusive wave system for an 
arbitrary perturbation. 

In the special case, when the scattering fidelity is measured from a reflection 
amplitude (a = b), we find 

/aa({/%}, K) 5 = N (t)~ l T a [*? (Z)l + 2 (Z P >/] , (24) 

where the geometric mean of the auto correlation functions, denoted by N(t) turns out 
to be time dependent. While the effective transmission coefficients are still zero and 
Z = 1, the auto correlation functions now read: 

/o»({«i}, M ; t) = faaiWjh {«;•} ; t) = ^(*) = T Q 2 [4t 2 (1) 7 + 2 (P )/ ] • (25) 

The integral (1)/ has been calculated in [25], with the result: At 2 = 1 — b 2 (t), where 
62 (t) is the two-point spectral form factor of the Gaussian orthogonal ensemble [23J. 
Hence, we obtain: 

f (1 1 X '1 A 4t 2 (Z) I + 2(ZP ) I 

This result will be used in section 13. 2\ where we discuss experimental results for 
perturbations due to the displacement of an S-wave scatterer. 
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3. Local vs. global perturbations 

A detailed discussion of the differences between local and global perturbations can be 
found in [22]. Consider fl2]), where a perturbation results in the change of several 
eigenvalues of the perturbation operator W. In order to affect the dynamics of the 
system (leading to the decay of the fidelity amplitude), one may either change only a 
few eigenvalues by a large amount (local perturbation) or very many eigenvalues by a 
small amount (global perturbation). Both cases are considered in this section. 

3.1. Global perturbation 

This was the first problem solved in the context of fidelity decay of quantum-chaotic 
systems [HI [9J [12] . Experimentally, the perturbation was realized in a chaotic microwave 
billiard by displacing one of the straight billiard boundaries. If described by (pQ) and (EJ), 
W a may represent absolute displacements with respect to some initial position. Then, 
its eigenvector representation 

N 

W^ = X>» \Vj)( v i\ ( 27 ) 

runs over a large number N of states. According to section 12. 3[ and in particular (ITS]) 
and ( fl6l) . the effective transmission coefficients become 

T _ 2tt iWj(P) -iwj(a) 

j ~ N [l + mvjjiP) /N] [1 - m Wj (a) /N] 

= 2m5 J (l-in5 J ) + 0( y [w J ((3)/N] 3 , [wj(a)/N] 3 ) , (28) 

where Sj = [vjj((3) —Wj(a)]/N. In this setting, global perturbations are characterized by 
the fact that the contribution of each individual term is negligible, while the perturbation 
becomes noticeable only because it is the sum of very many such contributions. This 
allows to perform a Taylor expansion of In Z with respect to the coupling parameters Sj. 
Starting from the Taylor expansion of Zj with respect to the transmission coefficients 

Zj = [1 - Tj (t - r)} [1 + 27} r + T 2 x 2 )^ 2 

= [1 - T 3 (t - r)] [1 - rT 3 - [x 2 - 3r 2 ) if/2] + 0(Tf) 

= l-tT j + [rt + r 2 /2 - x 2 /2] T 2 + <D(T*) , (29) 

we insert f)28p into the above Taylor expansion, and obtain 

\nZ = -2 7 ri^(5 J t-27r 2 ^(r 2 + (2r + l)t-t 2 -x 2 ) 5 2 

j j 

+0(K(/3)/iV] 3 , K-(a)/iV] 3 ) • (30) 

In order to obtain a well defined function Z(t, r, x) in the limit of N — > oo, and vanishing 
perturbation: 5j —> 0, the parameters 5j must scale with an appropriate negative power 
of iV: (i) For 5j ~ A^ _1 , InZ would converge to the finite value — 27ri <5j t. (ii) For 
Sj ~ iV -1 / 2 , the sum ^ . 5j could still converge, if the Sj had different signs. In addition, 



Scattering approach to fidelity decay 9 

the sum ^ . 5j 2 would always converge, while any higher order terms would vanish, (iii) 
For powers larger than —1/2, the sum £V 5 2 would always diverge, and the function 
Z(t,r,x) would not be well defined. Hence, the cases (i) and (ii) are the only viable 
options, where 

= E 5 / = A2 ( 31 ) 

3 3 

converge to finite values. This finally leads to 

lim Z = exp(-27ri 5 S t - tt 2 A 2 P ) , (32) 

N— >oo 

with P Q given in (IT7|) . Note that by taking the absolute value squared of the fidelity 
amplitude, the dependence on d~ s disappears and with it any possible problems with the 
convergence of this term. 

To conclude this section about global perturbations, let us discuss the random 
matrix model for fidelity decay, as it has been first introduced in [7]. This model may 
be written as 

H a = H + aV , (33) 

where the matrices Hq and V are independent GOE matrices, normalized in such a way 
that the mean level spacing in the centre of the spectrum of Hq is do = 1, while for the 
perturbation matrix it holds 

(Vij Vki) = 5jk5u + S ik Sji . (34) 

Now, representing H a in the eigenbasis of V, the perturbation becomes diagonal with 
eigenvalues Wj(a) showing a semi-circle distribution between —2ay/N and 2ay/N. Then, 
according to (j2J): 

«i = o, ^ = x w M =► s i = -*jT' ^ 

which is of order A r_1//2 , indeed. From the semi-circle distribution of the eigenvalues 
Wj(a) it follows that 

3 3 3 3 

This shows that (1321) applies for this case if we set 5 S = and A = a. From (123]) it then 
follows that 

which agrees precisely with the result obtained in [8] 



f ab (t) = (P e-^^) i , (37) 
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3.2. Local perturbations 

In the case of local perturbations, W a and Wp differ strongly in a relatively small 
subspace. A Taylor series expansion as in the previous case is therefore not useful, 
and it is also less likely that [W a , Wp] = 0. Thus, it seems necessary to redefine 
the perturbation by considering H a as the unperturbed system and Wp — W a as the 
perturbation. In doing so, it is assumed that including W a into H doesn't change 
its statistical (i.e. random matrix) properties. We may then choose a basis in which 
Wp — W a is diagonal, the transformation into that basis leaves the random matrix 
ensemble for the new H invariant, so that we arrive again at a description in which the 
perturbation is diagonal. 

In contrast to the previous case, we have now only a small number of non-zero 
diagonal elements, while each of them may be very large. In principle, each element 
alone can cause the fidelity to decay as fast as in the case of a global perturbation. In 
what follows, we restrict ourselves to the least biased case, where the states \a) and 
| b) coupled to the measurement channels [cf. (jSJ) and (jHJ)] are not involved into the 
perturbation. In this case, the results become again independent from the measurement 
channels, and one can then repeat the previous argument to show that scattering fidelity 
(if measured in transmission) and fidelity amplitude must coincide. 



3.2.1. Moving scatterer This case refers to the displacement of a small scatterer from 
a position f\ to another position f%. As explained above, this is modelled by the 
unperturbed Hamiltonian consisting of the system with scatterer at position fi, while 
the perturbation consists in removing the scatterer from position fi and placing it at 
position f 2 . For a point-like scatterer, the effect of the scatterer can be described by 
one single state (the perturbation operator corresponding to that scatterer has only one 
non-zero eigenvalue). Therefore, 

H a = H , Hp = H + w(P) ( \v 2 )(v 2 \ - h)H ) . (38) 

For the scattering approach to fidelity, this means that the perturbation must be 
described by two effective transmission coefficients: 

1 — m di 1 + i7r d\ JS 

Here, 8\ may become arbitrarily large so that T\ as a function of b~\ traces a path in the 
complex plane which starts at T\ — and ends at T\ — 2, while T 2 = T*. Then 

n r Tlit - r) - L7T (t ~J\ <«> 



inserted into fl23|) or fl26|) yields the exact analytical expression for the decay of the 
fidelity amplitude. 
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3.2.2. Comparison with experiment The perturbation we just described, applies 
precisely to an experiment published in [201 E]. There, a small disk of diameter 
4.6 mm has been moved in steps of | Ar| = 1 mm through a rectangular two-dimensional 
microwave billiard with 19 additional random scatterers. Then, the reflection spectrum 
has been measured for 300 different positions of the moving disk in a frequency range 
from 3.5 to 6 GHz. In this frequency range, it was still possible to extract resonance 
positions and amplitudes by Lorentzian fits. The statistical properties of the spectrum 
as well as the wave functions were in agreement with the random matrix expectation 
for quantum chaotic or weakly disordered systems. 

From Berry's model of the random superposition of plane waves [26], it is possible 
to obtain a connection between the displacement |Ar| of the movable disk and the 
parameter Si which measures the perturbation strength. Translating the corresponding 
equation from [21] to our system of units and parameters, we obtain 

*i = f V 7 ! " Mk\Ar\) 2 , (41) 

where a is a dimensionless factor related to the electromagnetic properties of the movable 
disk, which can be determined independently (e.g. from the variance of the level 
velocities). For the wavenumber k, we choose a value which corrresonds to the frequency 
in the centre of the range considered, which yields 
2tt f 

k = — - = 0.996 cm- 1 . (42) 
c 

For the displacements considered in [21], the relation between Si and |Ar| is still 
approximately linear, as can be seen from the fact that k |Ar| is small as compared to one 
in all cases (see the captions of figure [IJ) . Finally, we find that a = 1 provides the best 
agreement between the theory and experiment. Using time independent perturbation 
theory, the authors of [21] obtained for the decay of the fidelity amplitude the following 
asymptotic result: 

y/1 + (ASit) 2 

valid for finite Si t, in the case where Si — > and t — > oo, and in agreement with our 
asymptotic result (l6Tj) . derived below (section HJ). 

In figure [T] we show the experimental data for the decay of the absolute value 
squared \f aa (t)\ 2 as obtained in [2TJ. In the experiment, this quantity is obtained from 
ensemble averages of the respective correlation functions, introduced in (I14p . The results 
are compared to the theoretical predictions based on the perturbative result, fl4"3"j) . and 
on our exact analytical expression for a reflection measurement, (126]) . 

We focus here on the behaviour of the fidelity at small times, where the asymptotic 
formula is expected to be less accurate, and indeed, we find significant deviations for 
the cases Si ~ 0.14 (green points vs. dashed green line) and Si ~ 0.28 (blue points 
vs. dashed blue line). For these cases, the experimental fidelity decay has a notable 
linear component at small times, which cannot be reproduced by the perturbative 
formula, ( 1431) . By contrast, our exact analytical result contains that linear component 
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Figure 1. Experimental data for the fidelity decay due to a moving scatterer from j21j . 
compared to the approximate (dashed lines) and to the exact theory (solid lines). 
The different colours, red, green, and blue, correspond to different displacements 
Ar = 1 mm (Si ps 0.07), 2mm (Si w 0.14), and 4mm (Si ps 0.28), respectively. 



and therefore agrees much better with the experiment (solid lines). Still, some differences 
remain for 5\ ~ 0.28. We believe that these are due to problems on the experimental side. 
One error source consists in the wide frequency range used, which, according to (HIT) leads 
to a considerable variation in the perturbation strength. Another problem is related to 
the upper end of the frequency range, which implies rather small wavelengths, for which 
the scatterer to be moved may no longer be point like. For a more significant test of 
our analytical formula, one would need a different experimental design, providing higher 
accuracies at strong perturbations in the vicinity of the Heisenberg time. 

In the perturbative result, ( 1431) . it makes no difference whether the measurement 
is performed as a transmission measurement with two measurement antennas or as a 
reflection measurement with only one. As we have seen from (1231) and (|26|) . for the 
exact analytical result this is no longer true. In figure |2} we compare both cases for four 
different perturbation strengths. The results show that the difference \f aa (t)\ 2 — \fab(t)\ 2 
is clearly present, but usually quite small. For very large (5± = 0.56, narrow peak at 
t < to) as well as for very small perturbation strengths, it seems that the difference 
tends to disappear. This is consistent with our treatment of the perturbative case in 
Sec. IH 

4. Perturbative regime 

In the perturbative regime, the fidelity amplitude of ([!]) can be calculated using first 
order time- independent perturbation theory [2]. For {\j)} denoting the eigenbasis of Hp 
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Figure 2. Comparison between fidelity decay, measured in reflection \f aa (t)\ 2 and in 
transmission \f a b(t)\ 2 ■ Different colours correspond to different perturbation strengths. 
Si = 0.07 (red line), 0.14 (green line), 0.28 (blue line), and 0.56 (pink line). 



and V = H a — Hp, we then find 

fit) « ]T(a|j> e 2 *^)* e -^M+U\vmt {j \ a) = £ \(j\ a )\ 2 e-^Wtft ^ (44) 

3 3 

where \a) denotes some initial state, and {Ej(f3)} denote the eigenvalues of Hp. This 
expression shows that in the perturbative regime the fidelity decay depends on the 
product between time and perturbation strength. The perturbative result becomes 
exact only in the limit of vanishing perturbation strength. To yield a finite value for 
the fidelity amplitude, time must then tend to infinity such that the product between 
perturbation strength and time remains constant. We therefore define the perturbative 
regime as the limit 

t-Kx, Vj : Tj -> , (45) 

such that t ^2,Tj and t 2 Ylj T 2 remain both finite. In what follows, we calculate f a b(t) 
and f aa {t) in that limit, starting from the exact analytical expressions ( )23|) and ( 126|) . 
via an asymptotic expansion of the respective integrals. This is done in two steps. 

4-1. Step one 

Here, we will demonstrate that 

<■■■>'-<■■■>*. »" h <"^ = £ dr jf dB pfeijF^Py - («) 

where the ellipsis above may be replaced by either Z P Q or At 2 Z. Here and in what 
follows, the symbol ~ denotes the perturbative limit we are interested in. 
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For the first case, our claim follows from 

f* , f, (t-r)(r + l-t)P Z n 

/ dr du ; " ~ . (47) 

Jt-i A/f {2u+l){t 2 -r 2 + X 2 ) 2 V ) 

Since Pq > and < Z < 1 in the whole region of integration, it is sufficient to show 
that f )4T|) holds for Z — 1. Furthermore, since we can maximize (t — r)(r + 1 — t) in the 
interval t — 1 < r < t by 1/4, it is sufficient to prove that 

r , r 2 + (2r + l)t - t 2 - a; 2 

max / du — . . -r— ~ , (48) 

t-i<r<tj^ (2u + l)(t 2 - r 2 + x 2 ) 2 ' v ; 

where we have used that Pq = 2 (r 2 + (2r + l)t — t 2 — x 2 ). Denoting this integral with 
J, we realize that 

[ r r 2 + (2r + l)t - t 2 

J< J^ (2u + l)(t 2 -r 2 + x 2 ) 2 m 

Then, because r 2 + (2r + l)t - t 2 = t + 2r 2 - (r 2 - 2rt + t 2 ), 

J < f du , 2 ^ + f = (2t 2 + t) f du , 2U+ ,] . (50) 
7^ (2« + l)x 4 1 (2r + l) 2 « 4 1 ; 

Evaluating the last integral we finally obtain: 

t(2t+l) / %/T+l 3r + l\ 



(2r + l) 2 V 3r 3 / 2 3r 3 

which completes the proof. For the second case, we replace Pq with 2i 2 one arrives at 
the same result, which may be seen from (]50p . 

4-2. Step two 

According to (|46|) the perturbative limit only requires integration of u up to u = \fr. 
This implies that 

Z = T\ i-W-r) TT^^ (52) 

11 /l x9T.,± ^2 „ 2 11 A /l + 2T,-t 



since (t — r) is of order one, t — 1 < r < t, and 

T 2 x 2 = T 2 u 2 ^1 < T 2 r ^L+L „ o . (53) 
3 3 2m +1 J 2^+1 v ' 

Therefore, we obtain for f a b(t) = (Z Po)i- 

This simply follows from the fact that 1 = (Pq)i ~ (Pq)a- For the scattering fidelity in 
a reflection measurement, we obtain 
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since 62 (t) in (126 j) tends to zero for large times. Here, it only remains to treat the first 
term in the nominator: 

«^>-n 7r = 1 « 2 «-n 7 rTi!p- (56) 

since 1 = At 2 ~ At 2 (1)a- Thus, in the perturbative regime, we obtain the same 
result no matter whether we perform a transmission or a reflection measurement: 

faait) ~ Ut) ~ / pert (t) = J] . (57) 



^.5. Global perturbation 

Global perturbations are discussed in detail in Sec. I3.1[ where ( 128)) relates the effective 
transmission coefficients with the perturbation strengths 5j. Taking also (EJTj) into 
account, we may write for f(t) up to second order in the perturbation strength: 

In /pert(i) ~ - -^ln[l + 47ri^(l-i7r^)t] ~ - ^ [2n\5jt + 2n 2 5 2 {2t 2 + t) ] .(58) 

Since we are working in the perturbative regime, where t goes as fast to infinity as the 
Sj go to zero, terms containing 5 2 1 can be neglected. Finally, we obtain 

\nf pcTt (t)~ e - 2M -^ x2t2 . (59) 
4-4- Moving scatterer 

Inserting the effective transmission coefficients from ( 1391) describing a moving scatterer 
into (1ST)) , we find 

4i M 
1 

In the perturbative limit considered here, means that T\ — > 0, t — > 00 such that T\ t 
remains constant. This implies however that also b~\ — > such that 8\ t remains constant. 
Therefore 

/pert(t) - II -4iM~ x = ; 1 (61) 



/pert(0 



(60) 



5. Comparison with parametric level correlations 

For any quantum mechanical model of the form 

H(X) = H + A V , (62) 

we may consider the level dynamics obtained from plotting the eigenvalues of H(X) as 
functions of A. For convenience, we assume here again that for any value of A, the 
average level spacing is one. In a typical random matrix model, one would eventually 
choose V from a Gaussian random ensemble with non-diagonal elements of unit variance. 
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The parametric level correlations X(X,r) describe the probability to find two 
eigenvalues, one of H(0) and the other one of H(X) at a distance r. The quantity 
to be compared to the fidelity amplitude is the Fourier transform of the parametric 
level correlations [Tj 



K(X,t) = j dr e 2mrt [X(X,r) - 1] . (63) 

Note that for A — > 0, this quantity converges to the complement of the two-point 
form factor: K(0,t) = 1 — 62^) [23 |23]|ij] The relation discussed in [14] is a relation 
between the parametric level correlations on the one hand, and the fidelity amplitude 
f\(t) = f(X,t) on the other. It may be expressed as 

f ^ t) = ^dW) K{x,th (64) 

with (3 being the Dyson parameter [28] which is one in our case. We consider systems 
with an anti-unitary symmetry such as time reversal invariance. 

5.1. Global perturbation 
From [22] we find 

r r°° r 1 (x^Xi- u') 2 (i- u' 2 )e inr +( XiX2 -^e- n2x2p ° 

A- ( A, r) = 1+Be // dA ld A 2 1 v (AiA2 ( ; + ) 2 L;-Vai-,-; • < 65 > 

where 2P = 1 + 2Af A| — Xf — X 2 , — n' 2 . We will see below, that this quantity is 
precisely the same as Po defined in the previous section, in ( flTl) . The first substitution, 
ji' — > n = (A1A2 — //') /2, yields 



/■(AiA 2 +l)/2 4„ 2 n _ ,/2\ „2i7rr +At -7r 2 A 2 P 

X(A, r) = 1 + 2 Re / / dAidA 2 / d/x { - / / x2 x9 ^ . (66) 

'l 7(AiA 2 -l)/2 



(l + 2A 1 A 2 /i / -A2-A2-^: 



In order to shorten the expressions, we keep writing // which must be understood as 
being a function of \x. Now, we can switch to the Fourier transform, which turns the 
Fourier factors in delta functions: 



(a^+d/2 [8(t + fi) + 8(t-n)]4fi 2 (l 



/2\ „-7T 2 A 2 P 



This shows that the function K(X,t) is symmetric in time. In what follows, we 
thus assume t > 0. The remaining delta function already allows to eliminate the /1- 
integration. However, before actually doing so, we perform a variable transformation 
on the Ai, A 2 integrals: 

(Ai, A 2 ) -> (r', x') , r' = A X A 2 , x' = A 2 - Ai (68) 

I The definition of K(X,t) in [14] uses the wrong sign, while the final expressions for Ki(X,t) misses 
the variable v in the nominator of the integrand. 
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The Jacobian of this transformation is simply J = (Ai + A2) 1 = (V 2 + 4r') l l 2 . 
Therefore, 

A,i = / dr' / / d/i 

4/x 2 (1 - r' + 2/i) (1 + r' - 2/i) - //) e^ 2 p ° 

[1 + 2r'(r' - 2/i) - x' 2 - 2r' - (r' - 2/i) 2 ] 2 ^ ^ 

For the /i-integral not to yield zero, it must hold that (r' — l)/2 < t < (r' + l)/2. This 
modifies the limits of the r'-integral as follows: 

K(X t) = / 2,+1 dr'f '"-^ (l-^ + 2t)(l + ^-2t)e-»"> 

1 ' J J^v-th-r' V^TV [1 + 2r'(r' - 2() - ^ - 2r< - (r< - 2«) 2 P l ™' 

Further substitutions: r' = 2r + 1 and a;' = 2a; and the fact that the integrand is a 
symmetric function of x, yield 

K(\, t) =4? f dr rdx( t -r)(r + l- f )e-^^ aJ , ^ 

Comparing to (1371) . it is now easily checked that K(X, t) as defined here fulfils the fidelity 
amplitude - parametric form factor relation ( 1641) . 

5.2. General perturbation 

In [22] it is shown that parametric level correlations can be calculated for arbitrary 
perturbations. According to this reference, the term describing the global perturbation 
Ogiob = 7r 2 A 2 P must be replaced by a = cr glob + <t 1oc , where 

f 2i4A 1 A 2 - K ;. 2 (A 2 + A 2 - 1)" 



o-ioc(Ai, A 2 ,/i') = ^5Z ln 



i-i«>') 2 J ' 

j j 

Note that the discussion in Sec. 13.11 shows that the additional global perturbation 
could always be incorporated into a\ oc , via a large number of additional channels with 
infinitesimal perturbations. However, in order to establish the desired relation between 
fidelity decay and the parametric level correlations, it is important to have the parameter 
A describing the global perturbation at hand. 

Now, we should go through the calculation of K(X, t) again, replacing <7 g i b in (1661) 
with the more general expression a. As a consequence, the integrand is no longer 
real, which affects ( 167|) . While the delta function d~(t — fi) is multiplied with the same 
term as before, the second delta function 5(t + /i) is now multiplied with its complex 
conjugate. Therefore K(X,t) is no longer symmetric. Instead K(X,t) = K(X,—t)*, 
which nevertheless allows to continue the calculation without changes for t > 0. Only 
at flTTj) we need to express a\ oc in the current integration variables. That results in 

exp [ - <r loc (r, x,t)] = T\ 1 ~ T ^ t ~ r ) = , (73) 

2T jr + T?x 2 
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where Tj = 2i«'/(l + i /«'-), just as in (I2"8"j) . Inserting this expression into (17T|) and 
comparing to the general result fl23l) for the decay of the fidelity amplitude (a 7^ 6) we 
find that the following relation holds: 

/«*(«},*) = ^ *T({k;}A*)[_ o • (74) 

The fact that any global perturbation may always be modelled with a large number of 
additional terms in cri oc , allows for a final slight generalization: 



-1 







fabiWj}, A , t) = — — K({k'j}, A, t) 



4ir 2 1 2 

This relation constitutes our second important result. In practice, this relation means 
that one can obtain the fidelity amplitude, by measuring the change of the parametric 
level correlations under the increment of a global perturbation. 



6. Conclusions 



In this paper, we have used the results of [TH] to derive an exact analytical expression for 
the fidelity decay in a closed chaotic/diffusive wave system, under arbitrary Hermitian 
perturbations. For illustration, we used that result to re-derive the known formula for 
the fidelity decay in the case of a global perturbation [8J [9]. In a second application, 
we calculated the fidelity amplitude for a moving S-wave scatterer, and checked that 
it describes corresponding experimental results reported in [21] well. Finally, we 
generalized a relation between the fidelity amplitude and parametric level correlations 
introduced in [H] to arbitrary perturbations. 

In the present work, we restricted ourselves to matrix ensembles based on the 
Gaussian orthogonal ensemble (GOE). For the comparison with experimental data, this 
is the most important case. However, our results can also be translated to the Gaussian 
unitary ensemble (GUE), using the analog of the VWZ-formula published in [24]. For 
the Gaussian symplectic ensemble (GSE), the corresponding analytical expressions for 
the parametric level correlations and the correlations between scattering matrix elements 
are unfortunately not yet available, but we would still expect a similar relation to hold. 

It would be interesting to perform an experiment similar to the one analyzed 
in [201 EE], in order to verify our results with higher accuracy and for larger perturbation 
strengths. Particularly interesting would be the regime, where the perturbation strength 
depends in a non-linear way on the displacement of the scatterer, see (|4T|) . If the 
microwave experiment would allow to measure fidelity decay and parametric level 
correlations at the same time, one could test the applicability of the relation (1751) 
between both quantities practice. Finally, one may intend to generalize ( 1751) further 
to scattering systems and non-Hermitian perturbations (e.g. coupling fidelity). 
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